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Abstract 

We model the quantum confined Stark effect in AIN/GaN/AIN heterostructures grown on top of 
[0001]-oriented GaN nanowires. The pyro- and piezoelectric field are computed in a self-consistent 
approach, making no assumption about the pinning of the Fermi level, but including an explicit 
distribution of surface states which can act as a source or trap of carriers. We show that the pyro- 
and piezoelectric field bends the conduction and valence bands of GaN and A1N and transfers 
charges from the top surface of the nanowire to an electron gas below the heterostructure. As a 
consequence, the Fermi level is likely pinned near the valence band of A1N at the top surface. The 
electron gas and surface charges screen the electric field, thereby reducing the Stark effect. The 
efficient strain relaxation further weakens the piezoelectric polarization. We compute the electronic 
properties of the heterostructures with a sp 5 d 5 s* tight-binding model, and compare the theoretical 
predictions with the available experimental data. 
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I. INTRODUCTION 



Wide band gap nitride semiconductors are now widely used for light emission in the blue 
and ultraviolet range.- - - Thanks to large band offsets, GaN/AIN heterostructures are also 
promising candidates for fast telecommunication devices based on intersubband transitions 4 
or for high-temperature single photon emitters.- One of the specifics of nitride heterostruc- 
tures is the existence of large internal electric fields due to spontaneous polarization and 
strains (piezoelectricity).- These built-in fields might transfer charges in the devices, leading 
for example to the formation of two dimensional (2D) electron gases at the interfaces be- 
tween GaN and (Ga)AlN layers.- - — They might also separate the electrons from the holes in 
GaN quantum wells and Stranski-Krastanov (SK) quantum dots, thereby reducing the band 
gap and oscillator strength (quantum confined Stark effect).— - — 1 ^ It is therefore essential 
to understand and tailor the electric field in nitride heterostructures to suit a particular 
application. 

Whereas 2D layers usually feature a large density of dislocations, nitride nanowires 
offer the opportunity to make defect-free heterostructures thanks to the efficient strain 
relaxation.— ^ Single GaN quantum disks (QDs) between two A1N barriers have for exam- 
ple been grown on top of GaN pillars (20-50 nm diameter) with plasma-assisted molecular 
beam epitaxy (see Fig. [T]). 19 The exciton and biexciton luminescence of 1 nm thick GaN 
QDs has been observed, showing the potential of such heterostructures for nitride optoelec- 
tronic devices. 19 A strong red shift (below the bulk GaN band gap) has been subsequently 
observed for larger disk thickness, a signature of the quantum confined Stark effect.— The 
apparent electric field is however smaller than expected from a comparison with GaN/AIN 
quantum wells.— The effects of strain relaxation (decrease of the piezoelectricity) and finite 
lateral size on the electric field,— as well as the screening mechanisms are still unclear. 

Nitride SK dots have been modeled before with k ■ p or tight-binding approaches.— 1 ^ - — 
The key role played by the charges transferred by band bending, which screen the elec- 
tric field, has been emphasized in 2D GaN/GaAIN layers.- - — In this paper, we model the 
electronic properties of [0001]-oriented GaN/AIN nanowire heterostructures in an atomistic 
tight-binding framework.—^ We compute strains with a valence force field method,™ then 
the pyro- and piezoelectric field. We account for band bending with a semi-classical Debye- 
Hiickel approach, making no assumption about the pinning of the Fermi level, but including 
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a distribution of surface states which act as a source or trap of carriers. We show i) that 
the piezoelectric component of the field can be significantly reduced by the efficient strain 
relaxation in the nanowire geometry ; it) that the spontaneous and piezoelectric polariza- 
tions create an electron gas at the lower GaN/AIN interface and are likely large enough to 
create a hole gas in the upper A1N barrier. These electron and hole gases screen the electric 
field in the GaN QD and reduce the Stark shift ; m) that for carefully chosen dot and bar- 
rier thicknesses (realized experimentally) the GaN QD is empty at equilibrium, consistent 
with the observation of exciton and biexciton transitions. We discuss the magnitude of the 
electric field and the electronic stucture of the QDs as a function of the dimensions of the 
heterostructure. 

The paper is organized as follows: We review the methods in section [TTJ then discuss the 
electric field in GaN/AIN nanowire heterostructures in section HTT1 We introduce a simple 
ID model for the pyro- and piezoelectric field that reproduces the main trends. We finally 
discuss the electronic structure of the GaN QDs and compare our calculations with the 
available experimental data in section IIVI We analyze the dependence of the electronic and 
optical properties of the QDs on the geometry of the heterostructures. 

II. METHODOLOGY 

In this section, we introduce the methods used to compute the structural and electronic 
properties of the GaN/AIN heterostructures. 

Each nanowire is modeled as a 30 nm diameter and 150 nm long 30 cylindrical GaN 
pillar oriented along [0001], with the heterostructure on top (we assume metal-face polarity 
- the case of N-face polarity will be briefly discussed at the end of paragraph II Vp . The 
heterostructure consists of a lower A1N barrier with thickness t in f , a GaN quantum disk with 
thickness £qd and a upper A1N barrier with thickness t sup (see Fig. CE]). The dangling bonds 
at the surface of the nanowire are saturated with hydrogen atoms. 

The strains in the nanowire are computed with Keating's valence field model."' This 
model provides an atomistic description of the elasticity of tetrahedrally bonded semicon- 
ductors. It was originally designed for zinc-blende and "ideal" wurtzite materials with equal 
bond lengths and angles,— and has been recently adapted to arbitrary wurtzite materials 
such as GaN and A1N.— The elastic energy of the nanowire is minimized with respect to the 



3 



FIG. 1: (Color online) Structure of the GaN/AIN nanowire heterostructures. The radius of the 
nanowires is R = 15 nm. 

atomic positions using a conjugate gradients algorithm. The strains e a p on each atom are 
then calculated from the atomic positions using a method similar 33 to Ref. 
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The pyro- and piezoelectric polarization density is next computed from the strains on 
each cation (Ga or Al) : 

/ \ 



P = PnZ + 



2>&15 £ yz 



(1) 
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where z = [0001], P is the spontaneous polarization, and ei 3 , e 33 and e 15 are the piezo- 
electric constants of either GaN or A1N (see Tabled]). Poisson's equation for the pyro- and 
piezoelectric potential V p (r): 

Ko V r • K (r)V r V p (r) = V r ■ P(r) (2) 



is then solved on a finite difference grid^S (see appendix I of Ref. |27| for details), k is the 
dielectric constant (k — 9 inside the nanowire and k — 1 outside). 
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TABLE I: The material parameters for GaN and AIN^^-Sl 

The large pyre- and piezoelectric field E p = — W p in the heterostructure bends the 
conduction and valence bands and can therefore transfer charges from one part of the system 
to an other. It is for example known that the spontaneous polarization in GaAIN layers 
grown on GaN pulls out electrons from the GaAIN surface, which accumulate in a 2D 
electron gas at the GaN/GaAIN interface.- - — These electrons leave positive charges at the 
GaAIN surface, which can be ionized surface donors, emptied surface states, or even a hole 
gas. This redistribution of charges creates, in turn, an electric field opposite to E p , which 
can screen the latter to a large extent. 

The effects of band bending have been self-consistently computed in a semi-classical 
Debye-Hiickel approximation. The local density of electrons, n(r), and the local density of 
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holes, p(r) are calculated as:— 



n(r) = N c F 1/2 [-/3(E c -eV(r)-fi)} 
p(r) = iV v F 1/2 (£„ - eV(r) - //)] , 



(3a) 
(3b) 



where N c and are the effective conduction and valence band density of states of the 
material at point r, E c and E v are its conduction and valence band edge energies (see Table 
[I]), V(r) is the total electrostatic potential, and /i is the chemical potential or Fermi energy. 
Fu2 is the Fermi integral of order one-half and (3 = l/(kT), where T = 300 K is the 
temperature. We have, additionally, assumed that the nanowires were non-intentionally n- 
doped, with a concentration of donor impurities (silicon, oxygen or vacancies) Nd = 2 x 10 17 
cm~ 3 . The density of ionized impurities is:— 



where Ef, is the binding energy of the donor, which typically ranges from a few tens to a few 
hundreds of meV. 

As mentioned previously, surface states can play an important role in the electrostatics 
of nitride nanowires. They might act as a source^ - — or as a trap^ 1 ^ of carriers, effectively 
pinning the chemical potential in the band gap. Little is however known about the electronic 
structure of nitride surfaces.— On one hand, density functional theory (DFT) calculations 
on reconstructed GaN and A1N surfaces^ 3 - - — suggest the existence of occupied (donor-like) 
surface states above the valence band edge and empty (acceptor-like) surface states below 
the conduction band edge (as expected from simple considerations). On the other hand, 
the extensive literature about 2D electron gases in [0001] GaN/GaAIN heterostructures^ - — 
suggests the existence of dense (~ 10 13 cm~ 2 eV~ 1 ) surface donor states only ~ 1.5 eV 
below the conduction band of Gai-^ALN alloys (x ~ 0.4). Although the nature of these 
surface donors is still debated, oxygen has often been put forward.— It is not clear however 
that the same picture holds for non-polar Ga(Al)N surfaces and for surfaces of pure A1N, 
where the oxide is not the same. The situation is particularly tricky in nanowires, which 
expose different (polar and non-polar) surfaces. For the sake of simplicity, we assume in this 
work the existence of a uniform density of occupied surface states in the [E^E^] energy 
range above the valence band edge, and of a uniform density of empty surface states in the 
[E^jE^] energy range below the conduction band edge. The density of ionized occupied 
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surface states is therefore: 



1 + l e P[E+-eV(r)-n] 

N+ (r) = kT Df In 2 — T - — -j , (5) 

while the density of electrons trapped in the empty surface states is: 

1 4- 9p -P{Ei-eV(r)~n} 

N~(t) = kT DJ In - — -j . (6) 

D+ and D~ are the density of occupied and empty surface states, respectively (per unit 
surface and energy). The values of , E%, E± and E2 used in this work are also reported 
in Table HI They are representative of ab initio calculations,—"— and reproduce the pinning 
of the Fermi level on non-polar GaN surfaces.— We have varied = D~ between 5 x 10 12 
cm~ 2 eV _1 and 5 x 10 13 cm~ 2 eV _1 . Their effects will be discussed in paragraph IIHI We will 
show, in particular, that the electric field in the QD is weakly dependent on the model for 
the surface states up to large Df and D~ . 

In practice, the carrier densities n(r) and p(r) are computed on each Ga, Al and N atom, 
while the surface state densities iV s + (r) and N~(r) are computed on each hydrogen atom. The 
charge on each atom is then transferred to the finite difference mesh, and Poisson's equation 
for the total electrostatic potential V(r) is solved self-consistently with the Newton- Raphson 
method:^ 

/t V r ■ K(r)V r V(r) = V r ■ P(r) 
+ [ n (r)-p(r)-N+(r) + N;(r)-N+(r)]e. (7) 

The chemical potential /1 is adjusted to ensure overall charge neutrality of the nanowire. 

Finally, the electronic structure of the GaN QD in the potential V(r) is computed with 
a sp 3 d 5 s* tight-binding model.— - — In order to access the relevant states directly, a slice 
containing the GaN QD and 4 nm of each A1N barrier is cut from the nanowire. The bonds 
broken by this operation are saturated with hydrogen atoms, and a few conduction and 
valence band states are computed with a Jacobi-Davidson algorithm.—^ The convergence 
of the electronic structure of the QD with respect to the thickness of the A1N barriers has 
been checked. The above methodology has been implemented in an in-house code called 
TB_Sim. 
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t m { = 8 nm ; t QD = 4 nm ; t sup = 8 nm 
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FIG. 2: (Color online) The conduction band edge energy e c (r) = E c — eV(r) and the valence band 
edge energy £„(r) = E v — eV(r) along the axis of a nanowire with £qd = 4 nm and ii n f = t sup = 8 
nm. The reference of energy is the chemical potential /x = 0. The position of the A1N barriers is 
outlined in gray, and the top surface is at z = 0. 

III. THE ELECTRIC FIELD IN GAN/ALN NANOWIRE HETEROSTRUC- 
TURES 

In this section, we discuss the electric field in GaN/AIN nanowire heterostructures. We 
first analyze a particular case as an example. We then show that the electric field can be 
reproduced by a simple ID model in a wide range of dimensions. We finally discuss the 
main trends as a function of the geometry of the heterostructures. 

A. Example 

We focus as an illustration on a 30 nm diameter nanowire with a £qd = 4 nm thick GaN 
QD and t in{ = t sup = 8 nm thick A1N barriers. We set Df = D~ = 10 13 cm" 2 eV _1 . 

The conduction band edge energy e c (r) = E c — eV(r) and the valence band edge energy 
e v (r) = E v — eV^(r) are plotted along the axis of the nanowire in Fig. [3 The reference of 
energy for this plot is the chemical potential fi = 0. The top of the nanowire is located at 
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t m { = 8 nm ; t QD = 4 nm ; t sup = 8 nm 
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FIG. 3: (Color online) The bare (unscreened) and screened pyro- and piezoelectric potentials along 
the axis of the nanowire (£qd = 4 nm, ij n f = t sup = 8 nm). The position of the A1N barriers is 
outlined in gray, and the top surface is at z = 0. 

z = and the position of the A1N barriers is outlined in gray. 

The band discontinuities at the GaN/AIN interfaces are clearly visible. The heterostruc- 
ture undergoes a strong vertical electric field, which is almost homogeneous in the A1N 
barriers and in the GaN QD. The latter is empty of carriers (electrons and holes). The 
chemical potential is however pinned at the valence band edge at the top surface of the 
nanowire, and crosses the conduction band at the interface with the GaN pillar. Electrons 
therefore accumulate in the GaN pillar, while holes accumulate in the upper A1N barrier. 

This redistribution of charges follows from the pyro- and piezoelectric polarizations. 
Leaving aside piezoelectricity for the moment, the spontaneous polarization in GaN is 
P z = —0.034 C/m 2 , while the spontaneous polarization in A1N is P z = —0.090 C/m 2 . 
This polarization is equivalent to a distribution of charges o = —0.090 C/m 2 at the top 
surface, and a = ±(0.090 - 0.034) = ±0.056 C/m 2 at each GaN/AIN interface. Such a 
charge distribution, if unscreened, would create huge vertical electric fields and potentials 
of the order of 10 to 20 V in the nanowire (see Fig. [3]). 

The pyro- and piezoelectric field however bends the conduction and valence bands and 
tends to draw positive charges at the top of the nanowire, which screen the polarization. 

9 





FIG. 4: (Color online) The electrostatic potential V(r) in the (xz) plane containing the axis of 
the nanowire (£qd = 4 nm, i in f = t sup = 8 nm). The GaN and A1N layers are delimited by dotted 
lines. 

The potential actually rises the occupied surface states of the upper A1N barrier above the 
Fermi energy. They therefore empty, leaving positive charges at the surface and releasing 
electrons in the GaN pillar. At moderate electric field, the surface states would be able 
to provide enough charge to reach equilibrium, and the chemical potential would lie in the 
band gap at the top A1N surface. Here the electric field is however large enough to empty 
the N^t = Df(Ez — E±) = 2 x 10 13 cm -2 occupied surface states. The chemical potential 
then sinks into the valence band; a gas of holes forms at the top surface and provides the 
missing charges. 

According to this picture, the charge in the system is mostly distributed at the top A1N 
surface and at the GaN/AIN interfaces. As a consequence, the electric field is typical of 
a series of parallel plate capacitors, being almost homogeneous in the GaN QD and A1N 
barriers. This is further emphasized in Fig. HI which represents the electrostatic potential 
V(r) in a (xz) plane containing the axis of the nanowire. The equipotential lines are indeed 
parallel to the interfaces. The electrostatic corrections due to the finite cross section of the 
nanowire are therefore limited in the GaN QD and barriers in this range of dimensions. 

The effective density of states in the conduction and valence bands of GaN and A1N 
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FIG. 5: The ID model used for the analysis of the electric field in the GaN QD. 

are large enough to "lock" the potential at the interface with the GaN pillar and at the 
top surface (as a small variation of potential leads to exponential variations of the charge 
densities once the Fermi energy is in the bands). Hence, 

/i~E v (A\N)-eV(z = 0) (8a) 

at the top surface, and 

(i ~ £ c (GaN) - eV(z = -t hct ) (8b) 

at the interface z = —thct = — (^mf + ^qd + ^sup) with the GaN pillar. The voltage drop 
AV = V(z = —thct) — V(z = 0) across the heterostructure is therefore: 

eAV ~ E c (GaN) — E V (AYN) 



(GaN) + E„(GaN) - E V (A\N) 



(9) 



The voltage drop across the heterostructure is thus primarily defined by the band gap 
^ (GaN) of GaN and the valence band offset between GaN and A1N once the Fermi energy 
is pinned in the valence band of A1N at the top surface. The validity of this assumption will 
be discussed in the next paragraphs. 



B. A simple ID model 



We can derive a simple ID model for the electric field Eqd in the GaN QD from the 
above observations. For that purpose, we neglect finite size effects (R — >■ oo) and doping. 

11 



We assume that the polarization is homogeneous in the lower A1N barrier (P z = Pi n f), GaN 
quantum disk (P z = Pqd) and upper A1N barrier (P z = P sup , see Fig. [5]). This polarization 
is equivalent to a charge density <r sup = AP sup = P QD — P sup on the upper QD interface 
and (Ti n f = — AP in f = — (Pqd — Pinf) on the lower QD interface. We also assume that the 
difference of potential AV^ across the heterostructure is set by band structure effects [Eq. 
fl9]) or equivalent if other pinning of the Fermi level] . The electric field is then homogeneous 
in each layer and fulfill the continuity and integral equations: 

k k (P inf - P QD ) = AP inf (10a) 
«o« (P sup - Pqd) = AP sup (10b) 

tinf Pinf + ^QDpQD + t sup E sup = AV . (10c) 

We therefore get: 

£ QD = --^AP + ^ (11) 

where: 

^inf APinf ^supAP sup (12) 
tinf + ''sup 

is an average polarization discontinuity at the interfaces of the QD. Additionally, the electric 
field in the barriers is: 



Pinf = 

KqK 

E =J- 



AP inf - tinf + tsup AP 



ct 



AP sup - tinf + tsup AP 

thet 



+ — (13a) 

thet 
thet 



KqK 

The above equations hold as long as the QD is empty - which is also often desired 
experimentally. Neglecting quantum confinement in a first approximation, the QD is empty 
as long as the conduction band edge is above the Fermi energy, and the valence band edge 
below the Fermi energy throughout the dot. Assuming Pqd < (which is the case here), 
the QD is therefore free from holes if P„(GaN) — eV(z = — t sup — £qd) < /•*, and free from 
electrons if P c (GaN) — eV(z = — t sup ) > /i. Using Eqs. (JSJ), (EJ) and (|T3|) . these conditions 
respectively translate into the following constraints on t inf and t sup : 

P g (GaN) 

tinf < (14a) 

ePinf 

AV 

t SU p < . (14b) 
12 



Assuming fixed pyro- and piezoelectric polarizations, E- 1B £ and E sup are independent on t in f 
and t sup for given £qd and t^ et . Equations (fTl|) then show that the QD can be empty only in 
a finite range of positions within the heterostructure. The QD is indeed filled with electrons 
if it is too far from the surface, and filled with holes if it is too far from the pillar. Note, 
however, that quantum confinement will practically hinder the charging of the QDs by rising 
the electron and hole energies. The above constraints thus provide safe bounds for the design 
of nanowire heterostructures. 

Equations (TTTTl and (fl3"|) also give an estimate of the charge densities a s and a p accumu- 
lated at top surface and interface with the pillar, respectively. The continuity equation for 
the electric field indeed reads at this interface: 

«o« {E in f ~ E p n) = a p , (15) 

where E p n is the electric field in the pillar. Since E p n decreases rapidly away from the 
interface, 

Op ~ n nE ml . (16) 
Assuming that the tip of the nanowire is charge neutral at equilibrium, we then get: 

a s ~ - (tr p - AP inf + AP sup ) . (17) 

We can further split a s and a p into polarization and induced charges: 

a p = AP pil — n p e (18a) 
cr s = Pup + n s e , (18b) 

where AP p n = P pi \ — p n f, n p is the density of the electron gas at the interface with the pillar, 
and n s is the density of charges (ionized surface states+holes) at the top surface. The latter 
thus finally read: 

n p e ~ P P ii - P QD + tmf + tsup AP - k k^- (19a) 
n s e ~ -P QD + tinf + W AP - k k^ . (19b) 

£het ^het 

Note that t^et must be large enough for the electron gas to form at the interface with the 
pillar {n p > 0),- but this is usually not limiting the design of the heterostructure. 
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We can get a rough estimate of Eqd, n p and n s by neglecting piezoelectricity [P^ = 
Pup = Po(AlN), Ppii = P QD = Po(GaN)]. We then get from equations flTTJ and (fl9b|) : 

Pqd = ~ — tmf , + tsup [Po(GaN) - Po(AlN)] + ^ (20a) 

K K t-het ^het 



n p e = tini + tsup [Po(GaN) - P (A1N)] - k k— (20b) 

^het ^het 



t QD P (GaN) + (t inf + t aup ) Pq(MN) 

Tl s C — KqK . ( ZUC I 

^het ^het 

As a simple example, the pyroelectric field in a 4 nm thick QD embedded in an infinitely 
long nanowire (ihet — > oo) would be |Pqd| = 7.03 MV/cm. In a finite heterostructure 
with thct — 20 nm, the induced charges screen this field down to |Pqd| — 3.47 MV/cm. The 
density of the electron gas at the interface with the GaN pillar is then n p = 1.73 x 10 13 cm -2 , 
while the total density of charges (surface states+holes) at the top surface is n s = 3.85 x 10 13 
cm -2 . Therefore, the Fermi level is actually pinned in the valence band of A1N as long as 
the total (donor) surface states density is lower than iV+ it = 3.85 x 10 13 cm -2 . This critical 
density, although large, is yet not unreasonable for bare nanowire surfaces. We will however 
give further evidence in paragraph IIVI that the Fermi level is pinned at (or at least close to) 
the valence band edge of A1N. 

We have tested this simple ID model against the numerical solution of equations (J3])- 
(J7J). It gives excellent account of the electric field in the QD when ihet ~ 2P. The effects 
of the non-intentional doping are indeed negligible with respect to the amount of charges 
transferred by the pyro- and piezoelectric field. This model however tends to overestimate 
n p (as the electric field in the pillar actually decreases over tens of nanometers) and thus 
overestimates n s (by around 25% in the above example). Also, the piezoelectric polarization 
and field become inhomogeneous in thick heterostructures, as the strains are maximum at 
the interfaces and relax in between 2 ^ 1 ^ (see Fig. |6]). The ID model above is nonetheless 
very helpful in understanding trends and guiding the design of nanowire heterostructures. 



C. Discussion 



The amplitude of the electric field |Pqd|, computed with Eqs. (E])-(I7j) as the difference 
of potential along the QD axis divided by tQD, is plotted in Fig. [7K as a function of t m f 
and £ sup (£qd = 4 nm). As expected from Eq. (I20al) . the electric field increases with the 
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FIG. 6: (Color online) The conduction band edge energy e c (r) = E c — eV(r) and the valence band 
edge energy £„(r) = E v — eV(r) along the axis of a nanowire with tQ>r> = 4 nm and tjnf = t sup = 16 
nm. The reference of energy is the chemical potential \i = 0. The position of the A1N barriers is 
outlined in gray. The electric field is not constant in the barriers due to the inhomogeneous strains 
and piezoelectricity. 

total thickness t^et of the heterostructure, and ranges from ~ 3.5 MV/cm for thet — 12 nm 
to > 7 MV/cm for t^et = 36 nm. The electric field is slightly higher than expected from 
the spontaneous polarization, and does not fulfill the symmetry relation pQD(£inf, t 8np ) = 
EQ D (t sup , t inf ) due to piezoelectricity. This is further emphasized in Fig. [7)d, which represents 
the average AP obtained by inverting Eq. (TTTTl with the data of Fig. [7^. Three horizontal 
lines are also plotted on this figure for reference. AP pyro = 0.056 C/m 2 is the spontaneous 
polarization discontinuity at the GaN/AIN interface, which should provide a lower bound for 
AP. APcaN = 0.105 C/m 2 is the spontaneous and piezoelectric polarization discontinuity in 
a heterostructure biaxially strained onto GaN, and APain = 0.093 C/m 2 is the polarization 
discontinuity in a heterostructure biaxially strained onto A1N, which are the expected limits 
for thin and thick barriers, respectively. The actual AP lies between these bounds, as an 
evidence for piezoelectricity. The piezoelectric field, though still significant, is lower than in 
a 2D AIN/GaN/AIN quantum well, due to strain relaxation. The variations of AP result 
from a complex interplay between strain relaxation and charging (see discussion below). It 
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FIG. 7: (Color online) (a) The amplitude of the electric field |-Eqd| in the GaN QD as a function 
of and t sup (tQD = 4 nm). (b) The average polarization discontinuity AP deduced from (a) 
and Eq. (fTl]) . The two dots in parenthesis are charged with holes. 

is nonetheless worthwhile to note that a very good approximation to the electric field can 
be obtained with a constant AP ~ 0.077 C/m 2 (for given tgo and R) in a wide range of t m f 
and t sup . The value of AP slightly increases with decreasing £q D , up to AP ~ 0.082 C/m 2 
for t QD = 1 nm. 

As discussed above, the QDs might not be empty if they are are too far from the surface 
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or from the pillar. Eqs. (J3j) (JTJ) do not, however, properly take quantum confinement into 
account. We have therefore refined the assessment of the charge state of the QDs with 
the tight-binding model: we have tentatively assumed that the QDs were empty (setting 
N y = N c — in the dots so that they are free of carriers), and checked the position of the 
tight-binding band edges with respect to the Fermi energy. We find that all the QDs of Fig. 
[7] are actually empty, except those with t sup = 4 nm and t- m f > 8 nm, which are filled with 
holes. As expected from Eq. (I19bj) . the total charge density in the A1N barriers increases 
with t h ct, from n s = 2.39 x 10 13 cm" 2 for t in f = t sup = 4 nm, to n s = 3.54 x 10 13 cm -2 for 
^inf — t sup = 8 nm, and n s = 4.25 x 10 13 cm -2 for t in f = t sup = 16 nm. The Fermi energy 
is therefore pinned in the valence band of A1N at the top surface in all heterostructures 
considered here (n s > N^ t = 2 x 10 13 cm -2 ). 

We would finally like to discuss the role of the lateral surface states. The top surface 
states play a key role by releasing electrons in the GaN pillar, thereby screening the pyro- and 
piezo-electric field. The occupied lateral surface states of the upper and lower A1N barrier 
also act as a (secondary) source of electrons. Most of these extra electrons (as well as the 
donor electrons) are, however, trapped by the empty lateral surface states of the GaN pillar. 
As a consequence, the GaN pillar is effectively depleted far away from the heterostructure, 
and the Fermi level is pinned ~ 1.25 eV below the conduction band edge (see Fig. [2]).— This 
does not, however, have significant influence on the physics of the heterostructure. 



IV. ELECTRONIC PROPERTIES OF GAN/ALN NANOWIRE HETEROSTRUC- 
TURES 

We now discuss the electronic and optical properties of the GaN/AIN nanowire het- 
erostructures, and compare our results with experimental data. 

The tight-binding band gap energy E g of empty 4 nm thick GaN QDs is plotted as a 
function of the electric field Eq D in Fig. |HJ The corresponding values of t ini and t sup are 
reported between parenthesis. The excitonic correction is not included in this calculation 
and should further decrease the optical band gap by at most ~ 25 meV. The band gap 
energy is strongly red-shifted (below the bulk value) by the electric field (Stark effect). It 
depends almost linearly on Eqd and spans around 1 eV in the investigated range of ti n f and 
t sup . The ground-state electron and hole wave functions of a particular QD (t inf = t sup = 8 
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FIG. 8: (Color online) The band gap energy of GaN QDs as a function of the electric field (£qd = 4 
nm). The dotted line is a guide to the eye. The corresponding (t- m f, t sup ) are given (in nm) between 
parenthesis. The red square is the experimental structure discussed in the text. 
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FIG. 9: (Color online) The lowest hole (red/solid contour lines) and electron (blue/dashed contour 
lines) wave functions in a GaN QD {tqp = 4 nm ; t- m { = t sup = 8 nm). The gray dots are the Al 
atoms in the A1N barriers. 
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FIG. 10: (Color online) The calculated (TB) and experimental (Exp.) band gap energies of the 
GaN QDs. Experimental GaN/AIN quantum wells^ (QW) and Stranski-Krastanov (SK) quantum 
dots 16 with similar sizes are also reported for comparison. 

nm) are plotted in Fig. [9] As expected, the electron is confined at the upper interface, while 
the hole is confined at the lower interface of the QD. The electron and hole are, interestingly, 
both localized around the axis of the nanowire, where the strains are maximum, hence the 
piezoelectric field slightly larger than at the surface. This helps preventing one of the carriers 
from being trapped by the charged lateral surface states.— 
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Finally, we compare our theoretical predictions with the experimental results of Ref. 
In this work, the room-temperature luminescence of 1 to 4 nm thick GaN QDs embedded 
in 30 nm diameter nanowires showed clear evidence of the quantum confined Stark effect. 
We have therefore computed the electronic structure of a 1 nm thick QD (i in f = 10 nm, 



t S up = § nm), of a 2.5 nm thick QD (t 



inf 



9 nm, t sup 



10 nm) and of a 4 nm thick QD 



(iinf = 8 nm, t sup = 7 nm). The geometry and thickness of the barriers were chosen after a 
detailed analysis of the experimental TEM images. 56 The calculated and experimental band 
gap energies are plotted in Fig. [TUJ All the dots are empty, which is consistent with the 
observation of the biexciton in the 1 nm thick QDs. 19 The electric field ranges from 5.6 
MV/cm for tgrj = 4 nm to 7.3 MV/cm for tqrj = 1 nm. It is, as expected, much smaller 
than in GaN/AIN quantum wells^ (QWs) and Stranski-Krastanov^ dots (SKs) with similar 
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sizes due to to strain relaxation and screening by the electron gas and surface charges (see 
the comparison between QWs, SKs and QDs in Fig. [TOl) . 

The calculation reproduces the luminescence energies of the 1 and 2.5 nm thick QDs 
within the error bars, and the downward trend of the electric field with increasing QD size 
evidenced in the experiment. Still, the calculation underestimates the luminescence energy 
of the 4 nm thick QDs, where the Stark effect is the largest, by about 0.3 eV. Looking at Fig. 
[HI this suggests that the electric field in this QD is overestimated by ~ 20 %. The reason 
for this discrepancy is unclear at present. Increasing the density of surface states to pin 
the Fermi level in the bandgap of A1N increases the electric field and ultimately charges the 
dots. The experimental data have, moreover, been collected at low enough excitation power 
to prevent screening by photogenerated multiple electron-hole pairs.— i22,i£& The calculated 
electric field might be affected by the electromechanical coupling^^ (influence of the electric 
field on strains), by the uncertainties in the pyro- and piezoelectric constants of GaN and 
especially A1N, and by their dependence on strains (non-linear piezoelectricity).— 1 ^ A simple 
ID model however shows that the electromechanical coupling should reduce the electric field 
by at most ~ 5% (see Appendix[Aj. Although the nanowires of Ref. |20|are likely metal-polar, 
their polarity has not actually been assessed experimentally— Interestingly, we would like to 
point out that a N-face polarity would give a slightly better agreement between theory and 
experiment. In that case, the pyro- and piezo-electric field are reversed, so that a hole gas 
forms at the interface between the GaN pillar and heterostructure and electrons accumulate 
at the top A1N surface. The difference of potential across the heterostructure [Eq. (jHJ)] then 
becomes AV = E v (GaN) — E C (A\N) = —5.45 eV, which is slightly higher in magnitude than 
for the metal-face polarity [AV = E c (GaN) — E V (A\N) = 4.3 eV]. As a consequence, the 
pyro- and piezoelectric field in the QD are better screened, so that the luminescence in the 
4 nm thick QDs is raised by ~ 150 meV. Further experiments (for different barrier and QD 
thicknesses), as well as detailed polarity measurements might therefore be needed to get a 
complete picture. 

V. CONCLUSION 

We have modeled the quantum confined Stark effect in [0001]-oriented AIN/GaN/AIN 
nanowire heterostructures using a tight-binding approach. We have taken strain relaxation 
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and band bending into account in the calculation of the pyro- and piezoelectric field. We have 
shown that strain relaxation reduces the piezoelectric polarization, and that the electric field 
pulls out electrons from the occupied states of the top surface. These electrons accumulate 
in the GaN pillar below the heterostructure, thereby screening the pyro- and piezoelectric 
field. We suggest that the electric field is likely strong enough to pin the Fermi level in or 
close to the valence band edge of A1N at the top surface. As a consequence, the electric 
field is significantly reduced with respect to GaN/AIN quantum wells or Stranski-Krastanov 
quantum dots. This is in agreement with recent experimental data on GaN/AIN nanowire 
heterostructures.— The calculation however overestimates the electric field in thick quantum 
dots, which calls for further experiments with different geometries and detailed polarity 
measurements. We have, for this purpose, provided a simple ID model for the electric field 
to help the design of such heterostructures. We thank B. Daudin, B. Gayral, J. Renard and 
C. Bougerol for fruitful discussions. 



Appendix A: Electromechanical coupling 

In this appendix, we give an estimate of the electromechanical coupling and of its effect 
on the optical band gap of the nanowire heterostructures. 

We assume that the heterostructures are biaxially strained onto GaN. According to Fig. 
[7b, this should provide an upper bound to the electric field, hence to the electromechanical 
coupling. The in-plane strains e xx = e yy = £\\ in the GaN QD and A1N barrier are then: 

£|| = in GaN (Ala) 
en = 2.47% in A1N. (Alb) 



Following Ref. 



57, the vertical strain e 77 reads in each material: 



-2^ei, + —E z (A2) 
C33 c 33 



where C13 and C33 are the macroscopic elastic constants.— The first term is the "uncoupled" 
elastic strain while the second one describes the feedback of the electric field E z on the 
structure (electromechanical coupling). Solving Eq. flA2[) with Eqs. (|lip - (|13|) for the 
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electric field yields: 



Eqb = —8.15 MV/cm without electromechanical coupling (A3a) 
Eqd = —7.70 MV/cm with electromechanical coupling (A3b) 

in a 4 nm thick QD with t in f = 8 nm and t sup = 7 nm (experimental geometry). The elec- 
tromechanical coupling can therefore reduce the electric field by at most 5.75%. According 
to Fig. |H1 a 0.45 MV/cm decrease of the electric field would account for a ~ 120 meV 
increase of the bang gap. Additionally, the strains in the GaN layer are: 

e zz = without electromechanical coupling (A4a) 
e zz = —0.17% with electromechanical coupling. (A4b) 

Using the interband deformation potential a z = —11.3 eV in GaN,— the strains in the 
coupled system would further increase the band gap by ~ 19 meV. As a whole, the elec- 
tromechanical coupling can not, therefore, be expected to increase the band gap by more 
than ~ 140 meV. The same order of magnitude is obtained assuming the heterostructure is 
biaxially strained onto A1N. 
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